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ABSTRACT 

Results of a phenomenological model to estimate the GRB detection rate by the Fermi Gamma ray 
Space Telescope are reported. This estimate is based on the BATSE 4B GRB fluence distribution, 
the mean ratio of fluences measured at 100 MeV - 5 GeV with EGRET and at 20 keV - 2 MeV with 
BATSE, and the mean EGRET GRB spectrum for the 5 EGRET spark-chamber GRBs. For a 10% 
fluence ratio and a number spectral index ai — —2 at 100 MeV - 5 GeV energies, we estimate a rate 
of w 20 and 4 GRBs per yr in the Fermi Large Area Telescope field of view with at least 5 photons 
with energy E > 100 MeV and E > 1 GeV, respectively. We also estimate « 1.5 GRBs per yr in the 
Fermi FoV where at least 1 photon with energy £^ > 10 GeV is detected. For these parameters, we 
estimate w 1 - 2 GRBs per year detected with the Fermi telescope with more than 100 7 rays with 
E > 100 MeV. Comparison predictions for ai — —2.2, different fluence ratios, and the AGILE 7-ray 
satellite are made. Searches for different classes of GRBs using a diagram plotting 100 MeV - 10 GeV 
fluence vs. 20 keV ~ 20 MeV fluence is considered as a way to search for separate classes of GRBs 
and, specifically, spectral differences between the short-hard and long duration GRB classes, and for 
hard components in GRBs. 
Subject headings: gamma-rays: bursts — theory 



1. INTRODUCTION 

The Gamma ray Large Area Space Telescope, GLAST, 
was launched on 11 June 2008 and renamed the Fermi 
Gamma ray Space Telescope (hereafter Fermi) on August 
26, 2008. Fermi comprises two separate instruments, the 
Large Area Telescope (LAT), which is sensitive to 7 rays 
in the energy range from E fa 30 MeV to > 300 GeV, 
and the Fermi GLAST Burst Monitor (GBM), sensitive 
at 10 keV < < 30 MeV. At GeV energies, the Fermi 
LAT provides an increase in effective area over the En- 
ergetic Gamma Ray Experiment Telescope (EGRET) on 
the Compton Gamma Ray Observatory (CGRO) by al- 
most an order of magnitude, a smaller point spread func- 
tion, a field-of-view (FOV) larger by a factor w 5 than 
EGRET, and the capability to autonomously slew the 
LAT in response to a trigger from the GBM. The GBM 
uses Nal scintillators for triggering, like the Burst and 
Transient Source Experiment (BATSE), in addition to 
BGO (Bismuth Germanate) detectors for response in the 
- 0.2 MeV - 25 MeV energy range. The GBM views the 
unocculted sky, w 50% of the full sky, and the predicted 
long-duration GBM burst detection rate is w 150 - 225 
per year when triggering using standard BATSE detec- 
tion criteria (> 4.5(7 in at least two detectors in 1.024 
s in the 50 - 300 keV; iKippen et ali r20011. This can be 
compared with the BATSE detection rate of « 550 GRBs 
full-sky brighter than 0.3 cm'^ s'^ in the 50 ~ 300 keV 
band (Band 2002). 

Five EGRET 7-ray transients, coincident in time and 
direction with a BATSE GRB, were detected with the 
EGRET spark chamber at > 30 MeV wit h photon 
counts in excess of background (|Dinguslll995D . Knowl- 
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edge from the CGRO can be used to quantify expec- 
tations for Fermi. To simulate GRBs in the Fermi en- 
ergy band, we calculate the average fluence ratio between 
EGRET (100 MeV - 5 GeV) and BATSE (=i 20 keV - 
2 MeV) from the 5 BATSE GRBs that were detected 
with the EGRET spark chambe r. Using an averag e spec- 
trum consisting of an empirical iBand et al.l ()1993( ) spec- 
trum plus a power law, and the effective area of the 
Fermi LAT, we estimate the full-sky number of GRBs 
per year with different numbers of high-energy photons 
based on the BATSE 4B fluence distribution for long- 
duration GRBs (Paciesas et al. 1999). We make predic- 
tions for both long-duration and total GRBs; only the 
former class of GRBs were convincingly detected at 100 
MeV - 5 GeV energies with CGRO. Both Fermi and AG- 
ILE view 1/5*'' of the fuU sky at 100 MeV - 5 GeV 
energies, with the field-of-view defined by the solid angle 
where the effective area is > one-half the on-axis effective 
area. Thus we consider the realistic detection rate to be 
« 0.2 X the full-sky rate, without considering the poten- 
tial effects of autonomous repointings in the estimates. 

In § 2 we discuss the phenomenological model for GRBs 
based on BATSE and EGRET data, and estimate the 
uncertainty in GRB properties derived from EGRET 
and Fermi measurements. The predicted GRB detection 
rates and GRB contribution to the diffuse extragalactic 
background radiation are calculated in § 3. We note that 
these predictions were made in advance of the launch of 
Fermi (see arXiv:0807.0355 vl), so can be used to test 
our knowledge of GRB physics, in particular, whether 
additional hard, high-energy (> 30 MeV) 7-ray spectral 
components are present in GRBs. 

2. MODEL 

We develop a phenomenological model of GRBs based 
on the properties of EGRET spark-chamber GRBs, all 
detected with BATSE, and then make predictions for 
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the GRB detection rate with Fermi and AGILE. Based 
on BATSE and EGRET results, and from earlier GRB 
studies, the time-integrated (and resolved, where possi- 
ble) GRB photon number spectrum N(E-y) is assumed 
to be well-d escribed by the sum of the Band function 
(|Band et alJflQQa ). iVB(e), and a high -energy power law 
component, Ncie). We write this function using dimen- 
sionless notation e = Ej/nieC^ as 

Nie)^NBie) + NG{e) (1) 
where NB{e) takes the form 



+ k. 



e" exp[-e(a - /3)/e,J H{e;e 
, e"-'3exp(/3-a) H{e;e^„e 
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The parameters a and j3 are the low- and high-energy 
Band a and Band (3 indices, and i?br = WeC^Cbr ~ 100 
keV is the "break energy." The photon energy E'pk of 
the peak of the vF^ spectrum is related to E\y^ through 
the expression 



(2 + a)E^, 
a- 13 



(3) 



provided a > —2 and (3 < —2. The Heaviside function 
is defined such that H{x] y, z) — except at y < x < z, 
where H(x; y, z) = I . The term is the constant 
normalizing the number fluence to the > 20 keV BATSE 
energy fluence 20 keV) of a particular GRB, and 

is given by k^ = ^j^{> 20 keV)//^, where 



nieC^ 



(4) 



For model GRBs normalized to the BATSE fluence, 
we take the minimum energy ej^,^ = 30 keV, and maxi- 
mum energy e^^^ = 5 GeV. For the mean GRB spectral 
indices, we take a = —0.9, f3 = —2.2, and break en- 
ergy E^^^ = 250 keV, corresponding to the mean for these 
values from spectr al analyses of bright BATSE GRBs 
(jPreece et al.ll2OOO0 . 

For the high-energy spectral component, we assume 
that GRBs can be described by a power law at 7-ray 
energies written as 



Hie; e„,„,e_J 



(5) 



The average spectrum of GRBs measured by the EGRET 
spark chamber is ai = —1.95 ± 0.25 for 45 7 ra,ys de- 
tected with in 200 seconds of the trigger (jPingusI Il995l : 
iDingus et al. 1998), similar to the average EGRE T Total 
Abso rption Shower Counter (TASC) index (Catell Tet al.l 
I1998D . Band f3 values for BATSE GRBs have a large 
range, with most GRBs showing — 2 < /3 < 2.5, but lie 
in a different range than EGRET and Fermi emissions. 
The mean spectral index measured by EGRET for the 
five spark-chamber GRBs is not in good agreement with 
the index of BATS E and TASC GRBs at E > I MeV 
(jPreece et al.ll2000h . This suggests that the spectrum is 
probably not a simple power law at high energies, but 
different spectral shapes shoul d be considered (e.g ., cut- 
offs or extra components; see iKaneko et al] 120081 . and 
unusual spectral evo lution of high-energy components 
(jConzalez et al.ll2003D . 
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Fig. 1. — Approximate on-axis effective area of the Fermi LAT, 
EGRET, and AGILE. The symbols are the data, and the solid 
lines are the best-fitted functions. 



Since the spectral index distribution of GRBs at 
medium 7-ray energies is still debatable, which Fermi 
will clarify, in this paper we make predictions assuming 
that the > 100 MeV 7-ray spectrum can be described by 
a simple power law with ai = —2. For definiteness, we 
consider ai — —2 and ai — —2.2. The latter value corre- 
sponds to the softest spectrum within la of the EGRET 
measurements, and is the same as the Band /? used, effec- 
tively giving a single power law without break. Besides 
spectral index, our predictions are also strongly affected 
by EGRET uncertainties, which we will discuss in the 
next section. 

The term 

. _P$.(>20keV) 



I2 



where 



(7) 



is the normalization constant based on the average flu- 
ence ratio p between EGRET and BATSE. In our model, 
the uncertainties in the number of 7 rays that were 
detected by EGRET and the measured BATSE and 
EGRET fluences are absorbed into this factor p. Here 
e^jjjjj and e^^^ are the minimum and maximum energy 
at 100 MeV and 5 GeV, respectively; note that even 
though EGRET detected photons down to 30 MeV en- 
ergies, the effective area dropped rapidly below « 100 
MeV. The higher energy given by e^^^^ is set by the en- 
ergy where self-vetoing in EGRET reduces its effective 
area (see Fig. [T]). 

The number of s ource cou nts with energy e > e.y for 
different detectors (jPermer fc Dingus 2004 1 is given by 



/"OO 

S,{>e^)^fl A{e)N{e) 



de 



(8) 



where Ai{e) is the on-axis effective area of different de- 
tectors (e.g., EGRET, Fermi LAT, and AGILE) and 
is a collection function depending on the opening angle 
of the detector. Here we conservatively take the open- 
ing angle as energy-de pendent with = 68% contain- 
ment (jThompsonlll986l : lDermer fc Dingus! l2004[ ) though 
a larger region of interest (ROI) will not much change 
the results for bright GRBs with many 7 rays. The on- 
axis effective areas for EGRET (Ae), Fermi LAT (Ap), 
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TABLE 1 

Effective Area Parameters 



Effective 
Aroa 


aa 


ai 


12 


aa 


a4 


Ae 


-24.116 


36.435 


-18.000 


3.912 


-0.317 


Aa 


2.010 


0.501 


-0.080 






Af 


-11.583 


18.503 


-8.127 


1.555 


-0.109 



and AGILE {Aa) are take n from the best fit to the cah- 
brati on and analysis data (jPittori et al.ll20dl : iBelli et ahl 
120071 and references therein) as shown in Fig. [TJ The 
functions describing these areas are 



Ae{x) — 1^°'°^"'^^^°''^^' +"3^"' +"4^ 



(9) 



where x = logiJ, and parameters are given in Table 1. 
Most GRBs would take place at relatively large angles, 
~ 30 - 60°, from the Fermi axis, where the effective area 
is half this value. On the other hand, a larger ROI would 
include more counts. These represent two additional un- 
certainties in the rate estimates. 

The on-axis effective area of Fermi is shown in Fig. 
[TJ Signal analysis of photon events in the Fermi LAT 
depe nds on event r econstruction and background rejec- 
tion (jAtwood et al. f2009). The dijfuse class of events, for 
instance, minimizes background by rejecting tracks with- 
out good energy measurement in the calorimeter and re- 
quires both tracker and calibration information. GRBs 
can, however, be studied with the transient event class, 
which increases low-energy response at the expense of 
greater background using large regions of interest around 
the GRB direction. We consider only photons with 
E > 100 MeV, where the point spread function and con- 
sequently the background is smaller than for photons at 
lower energies, limited at higher energies to a maximum 
of 5 GeV. Over this energy range, both the EGRET and 
Fermi effective area is relatively flat, and measurement 
uncertainty can be more reliably estimated. 

3. FLUENCE ESTIMATES AND UNCERTAINTIES FOR 
EGRET GRBS 

In Table 2, we hst the pertinent data from EGRET to 
make estimates of Fermi detection rates. The first five 
columns of Table 2 give the GRB date-name, the angle 
of the GRB with respect to the z-axis of the EGRET 
telescope, the total number and the number of 7 rays 
with energies between 100 MeV and 5 GeV imaged by 
the EGRET spark chamber, and the maximum 7-ray en- 
ergy associated with the GRB. With this data, we can 
calculate the 100 MeV - 5 GeV fluence Fo.i_5 goV from 
an EGRET spark chamber GRB using the expression 



0.1-5 GcV — 



E 



A{E,,e) 



(10) 



with the sum being over the i = 1, . . . , iV^ 7 rays with 
energies from 100 MeV to 5 GeV detected by the EGRET 
spark chamber and associated with the GRB, located 
at the angle with respect to the pointing axis of the 
EGRET. 

By choosing a specific energy window between 100 
MeV and 5 GeV, we can minimize uncertainty in the 



GRB fluence estimate. This is where the EGRET re- 
sponse is best understood, as compared to the higher 
energ y regime with uncertain systematics (jStecker et alj 
l2008f ). or at lower energies where the point spread func- 
tion increases and the numbers of background photons 
are larger. Because the background photons are more 
likely to be at low energies, for example, due to difficul- 
ties in event reconstruction, the number of background 
photons, quoted below from the analysis of each specific 
EGRET GRB, is an upper limit to the number of back- 
ground photons in the 100 MeV - 5 GeV range. 

3.1. Fluence Error Estimate for EGRET Spark 
Chamber GRBs 

Several errors come into the estimate of the 100 MeV 
- 5 GeV burst fluence, including (i) the statistical error 
due to the finite number of detected photons; (ii) dead- 
time effects in the EGRET detector; and (iii) background 
spark-chamber events erroneously associated with the 
GRB. Other effects may also play a role, including (iv) 
errors in energy resolution (« 12 % and 9% 7-ray photon 
energy resolution errors at 100 MeV and 1 GeV, respec- 
tively), and (v) errors in effective area. For the EGRET 
effective area, we use the angle-dependent EGRET ef- 
fective area with the TASC in coincidence (Fig. 14 of 
Tho mpson et al.lll993f ). where the effective area for the 
associate burst (See Table[2]) is obtained by interpolating 
the data in Figure 14. 

For GRB 910503, 9 photons are detected within 15 s of 
the GRB trigger, and 6 are likely to be associated with 
this GRB, implying a possible background of 3 events 
(iSchneid et al.|[l99a). A 10 GeV 7 ray w as also detected 
84 s after this GRB (|Merck et al.lll995f) coincident with 
its location. Six of the 7 rays were detected within 2 
seconds, allowing as much as a factor of « 30% error for 
0.1 s EGRET deadtime if high-energy 7 rays are pro- 
duced roughly uniformly over this time interval. If such 
emission is produced in bursts shorter than the ~ 0.1 
s deadtime, then the GRB produce d a stronger fluenc e 
then measured. For GRB 910601 (Kwok et al.) I1993D . 
4 events were detected within 100 s of the GRB trig- 
ger, with 3 events greater 100 MeV 7 rays, when only 
1.5 events were expected. This represents a significant 
uncertainty on the fluence due to background. In con- 
trast, deadtime effects are unlikely to be important for 
GRB910601 assuming, as we do henceforth, that 7-rays 
are not produced in sub- 100 ms bursts. 

GRB 930131, the so-called "superbowl burst," is inter- 
esting in that it shows an i ntense ^0.1 s peak character- 
istic of a short hard GRB (jSommer "eTallfiggl . followed 
« 0.8 s later by a weaker « 0.4 s pulse, and also displays 
extended emissio n detected with BATE for 50 s after 
the GRB trigger (jKouveliotou et al.|[l99l ). Although its 
i9o(50 - 300 keV) = 19.200 ± 2.56 s, its t5n(50 - 300 keV) 
= 1.024 ± 0.091 s (jPaciesas et al.lll999D . so it is near the 
boundary between the lo ng soft and short hard GRBs 
(|Kouveliotou et al.l ll993). Note that short hard GRBs 
are o ften observed w ith Swift to produce extended emis- 
sion orris fc Gehr els 2008). Sixteen 7 rays are imaged 
within 25 s of the BATSE trigger, of which two pairs and 
a cluster of three arrive on a sufhciently short, < 0.2 s 
timescale that deadtime effects could have reduced the 
detected number of photons. Based on background ob- 
served before the trigger time of the GRB, only 0.04 
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TABLE 2 

EGRET Gamma-Ray Burst Spark Chamber Detections 



EGRET GRB 7 rays spectral index of —2 spectral index of —2.2 



Burst 
Date 


0a 




^^100 


Ema:c(Ge\) 


<T)<* 

BATSE 


<f>e 

EGRET 




EGRET 




<f>e 

EGRET 




910503 


24 


9» 


2 


10 


2.92{±0.01) 


1.28 


0.43 


2.68 


0.92 


2.20 


0.75 


910601 


12 


8 


3 


0.31 


1.5(±0.01) 


1.13 


0.75 


1.79 


1.19 


1.48 


0.99 


930131 


28 


18 


12 


1.2 


0.66{±0.11) 


26.6 


40.3 


24.8 


37.6 


20.3 


30.8 


940217 


10 


28 


10 


18 


6.6(±0.03) 


13.0 


1.97 


5.53 


0.84 


4.58 


0.69 


940301 


6 


7 


5 


0.16 


1.12(±0.01) 


1.14 


1.02 


2.48 


2.21 


2.05 


1.83 



Note. — "^Anglc from axis in degrees 
''Number of 7 rays imaged by EGRET spark ciiamber 

^Number of 7 rays imaged by EGRET spark chamber between 1 00 MeV and 5 GeV 
''BATSE 20 keV - 2 MeV fiuence, units of 10"* org cm"^, from lMeegan et al.l | |1996[) 



f r. - 



'I'EGRET /*J'B ATSE 



3 Nine events within 15 s after GRB trigger l lSchneid et al.lll992l) 

event s are expected by chance within 25 s (jSommer et aLl 
Il994[ ). so background should not be important for this 
GRB. This GRB was very intense at BATSE energies, 
but because most of this emission appeared in a first 
very short, very intense pulse, the BATSE 20 keV - 2 
MeV fluence is smallest of th e 5 sp ark-chamber GRBs. 

GRB 940217 (|Hurlev et al.lll994 . with spark-chamber 
emission extending to > 100 min after the GRB trigger, 
has associated with it the highest energy photon of all 
EGRET GRBs, w 18 GeV. A total of 28 7-ray photons 
were imaged by the EGRET spark chamber from this 
GRB, with 10, 8, and 10 appearing during the w 90 
s episode of strong BATSE emission, after the BATSE 
emission and before Earth occultation occurring at « 900 
s after the GRB trigger, and after the GRB appears from 
Earth occuhation at ~ 3700 to 5400 s after the GRB. 
For these three periods, respectively, 0.39, 1.8, and 2.9 
background 7 rays are expected. Because the GRB was 
so extended in time and the 7 rays arrive one or more 
seconds apart, deadtime effects can be considered to be 
small for this GRB. The duration of Earth occultation, 
representing > 62% of the known period of activity of 
GRB 940217, therefore represents a major uncertainty 
in the flucncc measurement of this GRB. 

Seven photons were i maged within 21 s of the GRB 
940301 BATSE trigger (jSchneid et al.lll995l ). though ap- 
pearing sparsely distributed in time so that deadtime 
effects are probably negligible. Background is also likely 
to be small for this GRB, as no other 7 rays are detected 
over a period of 150 s around the time of the GRB except 
for the 7 photons mentioned. 

For completeness, two other GRBs from the CGRO 
era should be mentioned. This includes GRB 910814, 
with 7 rays observed up to w 50 MeV with the TASC 
(when the spark chamber was disabled due to Earth oc- 
cultation; in any c ase this event was at 38.7° from the 
EGRET boresight: iKwok et al.lfT99l . The > 20 keV 
fluence of GRB 910 814 was ^ 2.51 x lO"'^ ergs cm^^ 
(jMeegan et al.lll996[) . Late in the CGRO mission, when 
the EGRET spark chamber gas was mostly depleted, 
GRB 990123 was observed with COMPTEL to « 18 
MeV. This was a bright GRB with a > 20 keV flu- 
ence = 5.1 X 10 ""^ ergs cm~^ (|Briggs et allll999| ). At 
redshift z = 1.6 (jKulkarni et al.lfr999( l. it provided the 
first confident lower limits on bul k Lorentz fac tor us- 
ing internal 77 opacity arguments (Baring '2006). Joint 
analysis of BATSE and EGRET-TASC GRBs indicates 
excess, > 1 MeV emission, though not necessarily ex- 
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BATSE 4B Catalog 
long-duration bursts 



Log [Fluence] (ergs cm ^) 

2. — Histo gram of BATSE 4B c atalog GRB fluence distri- 
data from lPaciesas et al.l I I1999I1 . 



tending to > 100 Me V, in bright, long duration GRBs 
(jGonzEilez et al.ll2009[ ). The > 20 keV fluence of the most 
notable of these bursts, GRB 941017, was = 3x 10""* ergs 
cm~^. 

Comparison of the distribution of the flucnces of GRBs 
shown in Fig. [2] in the BATSE 4B catalog with the flu- 
ences of the EGRET spark chamber GRBs in Table 2, 
and notable BATSE- TASC GRBs in the preceding para- 
graph, shows that the high-energy emission is, for the 
long duration GRBs, always associated with the most flu- 
ent GRBs. This is characteristic of a sensitivity limited 
detection of high-energy 7-ray emission. Consequently, 
we can use the increased capabilities of the Fermi tele- 
scope with respect to EGRET to make straightforward 
predictions by scaling from the EGRET results. Note 
that the weakest GRB in Table 2 in terms of fluence — 
GRB 930131 — is also possibly a short, hard GRB. 

3.2. Error from Finite Photon Number 

We consider a major error on fluence measurement to 
result from statistical fluctuations due to the small num- 
ber of photons imaged by the EGRET spark chamber and 
associated with the GRB, particularly since we limited 
the spark chamber events to those with energy between 
100 MeV and 5 GeV. It is a simple matter to make this 
error estimate using a Monte Carlo simulation, and then 
add other errors to get a complete error estimate. 

Let the integrated GRB photon number spectrum be 
denoted by dN/dEdA. The fluence measured between 
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energies Ei and E2 is given by 

F{Ei,E2) = / dE E 



\dEdA J 



The number of source photons is given by 
Ns{Ei,E2) = / dEA{E, 



El 



\dEdA 1 



(11) 



(12) 



The probability of detecting a photon of energy E is sim- 
ply given by inverting the expression 



J^^ dE A{E, 0){dN/dEdA) 
J^' dE A{E, 0){dN/dEdA) 



(13) 



for E, where r is a random number uniformly distributed 
between and 1. 

In the EGRET (and Fermi) energy range, A{E, 9) = 
const within a factor of w 1.8 for 100 < £;(MeV) < 5000 
and 6 < 30°, and can be approximated by a power law 
over this range, so we write A{E,d) = AqE~^.^ The ef- 
fect of changing 6 on A{E, 9) over the duration of the 
GRB is negligible for EGRET, which was a pointing 
telescope, but can be important in Fermi observations 
of GRBs. More accurate numerical calculations can be 
made following the analytic inversion given here. 

Assuming that dN/dEdA oc E^"'^ , the inversion of eq. 
(fT3|) is trivial, and we obtain for the detected energy the 
value 



(14) 



E{r) = [El-^ + r{El-^ E^-^] '^^'^^ , 

where s = u + ai. The mean photon energy in the energy 
range £'1 < i? < £'2 is 



{E) = 



E. 



2-ai 



E: 



1 — ai 



E. 



1 — ai 



1 



"1 



(15) 



The error in fluence due to the finite number, N^, of 
photons detected is then given by the root mean square 
deviation of the detected photon energies from the aver- 
age photon energy, that is. 



T.f^AE{n)-{E)f 



AE 



(16) 



The fractional fluence error, AF/F, as a function of 
the number of detected photons A'^^(100 MeV - 5 GeV) 
in the photon energy range between Ei — 100 MeV to 
E2 — 5 GeV range is shown in Fig. [3] for the simplest 
approximation of a flat effective area, so that u = 0. 
Because the fluence is weighted by energy, and we have 
restricted the photon energy to a finite range, a spectrum 
with ai = 2 gives the largest fluence uncertainty, which 
is well described by the function 



AF 



1.45 



(17) 



where that the factor of 1.45 is derived from fitting the 
data (with p = 2) in Figure [H We can take this ex- 
pression as providing the maximum error in fluence from 

^ See Fig.[T]and www-glast.slac.stanford.edu/software/IS/glast_ 
lat_performance.htm; restricting the photon energies to > 200 MeV 
makes this an excellent approximation. 
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Fig. 3. — Fractional uncertainty in measurement of the 100 MeV 
— 5 GeV 7-ray fluence due to finite number of photons from the 
GRB, for an underlying photon spectra described as a power-law 



with index ai = p. Calculation uses N ■- 
chosen randomly, except where noted. 



1000 sets of N.y photons 



the statistical uncertainty due to the finite number of 
photons. If there is a priori knowledge that the intrinsic 
spectrum has a specific number index, e.g., ai = 2.2, 
then the numerical coefficient of eq. (17) becomes rather 
smaller, « 1.4. Eq. (17) is a conservative estimate of 
fluence uncertainty, though it still omits other possible 
systematic uncertainties. Note that in the EGRET era, 
GRB 930131 had the largest number of photons, A^^(100 
MeV - 5 GeV)« 12, 7 rays, in the 0.1 - 5 GeV range. 
As is apparent from Fig. [31 finite photon number gives 
the major source of error to the fluence measurement of 
EGRET spark-chamber GRB. 

In comparison, the error associated with relative en- 
ergy uncertainty is, as already noted, at the ~ 10% level 
for EGRET. For the Fermi Gamma ray Space Telescope, 
the relative energy uncertainties at photon energies of 
0.1, 1, and 10 GeV are « 17%, 9%, and 8% for nor- 
mal incidence, and 15%, 9%, and 5% at 60° off-axis 
(|Atwood et al.ll2009l ). Because this can be a calibration 
uncertainty affecting all photon energy measurements 
systematically rather than statistically, this uncertainty 
on the fluence will start to dominate for those rare GRBs 
detected with Fermi with more than ~ 50 photons be- 
tween 100 MeV and 5 GeV. Effective area uncertainties 
also appear at the ~ 10% level. Consequently, the un- 
certainty from the flnite number of photons dominates 
fluence measurement for EGRET GRBs. 

From the discussion of EGRET GRBs in Section 3.1 
and their properties in Tabled we can now reliably esti- 
mate the fluence uncertainty associated with finite num- 
ber of photons, background and deadtime effects mea- 
sured with EGRET in the 0.1 - 5 GeV energy range. 
For GRB 910503, we take A^^(100 MeV - 5 GeV) = 2, 
essentially giving 100% uncertainty in the measurement 
of fluence (taking twice this value would give a 2a upper 
limit). For GRB 910601, we again conservatively take 
N^{100 MeV - 5 GeV) = 2 in view of possible back- 
ground. For GRB 930131, we take A^^(100 MeV - 5 
GeV) = 12, keeping in mind that a deadtime correction 
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Fig. 4. — Diagram plotting fluence measured between 100 MeV 
and 5 GeV with the EGRET spark chamber vs the 20 keV - 2 MeV 
fluence measured with BATSE for the 5 EGRET spark chamber 
GRBs. The EGRET fluence for GRB 940217 could easily be a 
factor of two larger due to Earth occultation during the active 
phase of this GRB. 

in the absolute fluence by tens of percent or more could 
be required, depending on the distribution of 7-ray ar- 
rival times, which Fermi will measure. For GRB 940217, 
we take 7Vt,(100 MeV - 5 GeV) = 8, subtracting two 
background photons that could have contributed during 
the first prompt and early afterglow in this GRB (note 
that no EGRET photons in the 100 MeV - 5 GeV range 
were detected after Earth occultation). For GRB 940301, 
we take 7V^(100 MeV - 5 GeV) = 5. 

3.3. Fluence Estimate from EGRET Photon Event List 

The statistical errors associated with the finite number 
of detected photons, being the largest source of error for 
the EGRET GRBs, are used with the values in Table [2] 
and eq. (fTO|l to produce Fig. [H This figure displays the 
100 MeV - 5 GeV fluence measured with the EGRET 
spark chamber vs. the 20 keV - 2 MeV fluence measured 
with BATSE. Bolometric corrections to a larger energy 
range are usually small if -Bpk, the peak photon energy 
of the vFi, spectrum, falls in this waveband. This is 
not always the case for these GRBs, and comparison of 
this figure with reported results in the Fermi era must 
consider bolometric corrections. 

This fluence-fluence diagram provides, potentially, one 
way to discriminate between different classes of GRBs on 
the basis of their 7-ray radiative efficiency. A much more 
accurate comparison will be provided for those GRBs 
with known redshift. Because deadtime and background 
effects have been corrected from these EGRET spark 
chamber GRBs, the fluence-fluence diagram gives a pic- 
ture where the long-duration GRBs have a radiative ef- 
ficiency p in > 100 MeV 7 rays of « 1% or a few % of 
the bolometric energy output. Assuming no class bias 
in the values in Table [21 these give the average deduced 
radiative efhciency (p) = 9% (column 8). This suggests 
that we use as a mean fluence ratio a number between 
« 10, before considering broad ranges of value of p or 
even separate classes with different p values. 



4. ALTERNATE FLUENCE ESTIMATE 

As a check on p, we also make an independent esti- 
mate using the number N^{> 100 MeV) of > 1 00 MeV 
7-raY photons that were detected by EGRET (jPingusI 
|1995() , assuming a hard component of the spectrum with 
spectral indices ai = -2 and ai = -2.2. The EGRET 
fluence in the energy range of 100 MeV to 5 GeV can be 
expressed as 



EGRET 



A^^(> 100 MeV) 



rSGoV pl-ai ip 
JlOOMcV^ "-^ 



(18) 

where A{E, 9) is the angle-dependent EGRET effective 
area obtained by interpolating the data in Figure 14 of 
[Thompson et alj ([1993,) . From Tabled we flnd that the 
average measured fluence ratios p for the 5 spark chamber 
bursts are p « 9% (column 10) and p ~ 7% (column 12) 
for spectra indices of ai — —2 and —2.2, respectively. 
This is consistent with our previous estimate of « 9%, 
because that referred to the > 100 MeV-5 GeV fluence 
to 20 keV-2 MeV fluence ratio, rather than the > 30 
MeV-5 GeV fluence to 20 keV-2 MeV fluence ratio, the 
latter of which is expected to be ~ 50% larger. However, 
these average fluence values of 9% and 7% are dominated 
by the boundary of long-soft and short-hard GRB 930191 
with the fluence ratios of 40% and 30% for spectra indices 
of -2 and -2.2, respectively. Four other bursts are long 
GRBs with the average fluence ratios of 1.3% and 1%. 
Hence, we expect the average fluence for the long GRBs 
to be larger, conservatively, we assume values of 5% and 
4% for spectra indices of ai = —2 and —2.2, respectively, 
to represent long-duration GRBs. 

However, as previously noted, the total 7-ray flux is 
probably larger than EG RET measured. The 7-ray emis- 
sion from GRB 940217 (|Hurlev et al.l |1994[ ) likely per- 
sisted during the period of Earth occultation, easily in- 
creasing the number of phot ons by ^ 2. Deadt ime ef- 
fects, as seen in GRB 930131 (ISommer et al.lll994D . could 
have caused the number of 7 rays to be undercounted. 
Besides, there is the examp le of GRBs with anomalous 
hard tail s, like GRB 990104 dWren et al.l l2dol) and GRB 
941017 ([Gonzalez et al.ll2003[ ). the latter with fluence 
in the 1 - 100 MeV range comparable to or greater 
than the sub-MeV fluence. Further, there is the low- 
significance M ilagrito detection o f the X-ray flash-type 
GRB 970417a (jAtkins et al.ll2003f ). which would not only 
have to be at low redshift z < 0.3 but also have much 
greater TeV than MeV fluence. Hence we expect the av- 
erage fluence ratios p between EGRET and BATSE to 
be « 10% and « 8% for spectra indices of ai = —2 and 
ai = 2.2, respectively. 

We use p = 0.1 for model with ai — —2 and p = 0.08 
for model with ai ~ —2.2 to define the value of ka. 
The estimated total number of long-duration GRBs that 
get detected by any det ectors are based on the BATSE 
4B fiuence distribution ([Paciesas et al.|[l999[ ). shown in 
Fig. [21 assuming an annual GRB rate of 670 GRB/yr, 
including untrigger ed GRBs wi th photon fluxes smaller 
than 0.3 cm-^ s'^ ([Bandl[2002[) . The BATSE 4B Cata- 
log contains 1292 bursts total, including short and long 
duration GRBs. There are a total of 872 long-duration 
GRBs that arc idcntiflable in the BATSE 4B Catalog. 
When we estimate the number of source counts over one 
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TABLE 3 

Fermi and AGILE GRBs/yr FOV (« 1/5*'' ) Predictions 



spectral index of -2 spectral index of -2.2 
Threshold Energy Fermi AGILE Fermi AGILE 

> 30 MeV (> 5 photons) 30" (18)'' 7" (4)'' 12" (O)"* 2" (< 1)'' 

> 100 MeV (> 5 photons) 20 (12) 2 (< 1) 8 (4) < 1 (< 1) 

> 1 GeV (> 5 photons) 4 (< 2) < 1 (0) < 1 (0) < 1 (0) 

Note. — This column has value with flucncc ratio of "10%; ''5%; "8%; "^4%. 




1 5 10 50 100 500 1000 1 5 10 50 100 500 1000 

Number of photons detected Number of photons detected 




1 5 10 50 100 500 1000 1 5 10 50 100 500 1000 

Number of photons detected Number of photons detected 



Fig. 5. — Prediction for the number of long-duration GRBs per year full-sky Fermi LAT (a and b) and AGILE (c and d) will see based 
on a measured fluence ratio extrapolated from a —2 photon spectum from EGRET/BATSE to EGRET LAT energies. The curves describe 
the burst rate assuming a —2 photon spectra for different energy thresholds. 



year period, the ratio 872/1292 is applied to represent 
the long-duration GRB predictions. 

5. RESULTS 

In Figs. [5] and [6l we plot full-sky Fermi LAT and 
AGILE GRB detection rates, assuming that the Band- 
spectrum extends to > 100 MeV energies. This simple 
model predicts that « 30, 20, and 4 GRBs will be ob- 
served per year with the Fermi LAT with at least 5 pho- 
tons with energies > 30 MeV, E > 100 MeV, and 
E > 1 GeV, respectively. This is assuming a hard spec- 
tral index ai — —2 and a fluence ratio p = 10%. For 
ai = —2.2, the model predicts that w 12, k, 8, and prob- 
ably no GRBs per year should be detected with at least 
5 photons with E > MeV, E > 100 MeV, and E > I 
GeV, respectively, now assuming that p = 8%. 

The analysis also predicts that k, 7 and « 2 GRBs 
per year with at least 5 photon counts at > 30 MeV 
and E > 100 MeV, respectively, should be detected by 



AGILE assuming ai = —2 and p = 10% (see Fig. 5d). 
For ai — —2.2, the model predicts that w 2 and proba- 
bly no GRBs with at least 5 photons with > 30 MeV 
and E > 100 MeV, respectively, will be detected by AG- 
ILE per year with p — 8% (see Fig. 6d). A consistently 
smaller number of GRBs/yr full-sky is predicted for both 
Fermi and AGILE for the steeper spectral indices, as ex- 
pected. In Table [3l we summarize the predicted number 
of GRBs per yr with at least 5 photons greater than a 
specific 7-ray photon energy in the FoV, « 2.5 sr, of the 
Fermi LAT and AGILE, considering GRBs with hard 
spectral indices ai = —2 and ai = —2.2, and for fluence 
ratios p = 10%, 8%, 5%, and 4%. 

These models assume that all GRBs have power law 
with ai w —2 and 10% fluence ratio or —2.2 and 8% 
fluence ratios, and these assumptions would be first to be 
examined if Fermi LAT shows strong violations of these 
predictions. For ai ^ —2 and p = 10%, we estimate from 
Fig. 5 that « 1 - 2 GRBs per year will be detected with 
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Fig. 6. — Prediction for the total number of GRBs/yr full-sky that 
spectra as described in Fig. 3 for an assuming —2.2 photon spectra. 

the Fermi telescope that have more than 100 7 rays with 
E > 100 MeV. The models also predicts that the GBM 
fluences of long-duration GRBs detected with the LAT 
should be large, > 10~^ ergs cm~^. The distribution 
of fluences can be compared wi th predicted distribu tions 
for Swift and the Fermi GBM (|Le fc Dermerj [20071) . 

From this treatment of the fluence ratio, we can es- 
timate the diffuse extragalactic 7-ray background pro- 
duced by GRBs using the 4B BATSE fluence catalog 
based on the EGRET/BATSE fluence ratio. From the 
4B Catalog there are a total of 1292 GRBs, of which 
872 bursts are long-duration GRBs. The total average 
fluence for the long-duration GRB is about « 1.40 x 
10~^ ergs cm~^. Assuming 670 GRBs/yr, we estimate 
that less than 1% of the diffuse extragalactic 7-ray back- 
ground could come from GRBs assuming « 8% or 10% 
fluence ratio. Unless there is a large class of p « 1 GRBs, 
which were not found with the EGRET spark chamber, 
this suggests that GRBs gives very little contrib ution to 
the diffuse extragalactic 7- ray background (e.g.. iDermeil 
[200I iCasanova et al.l[2007^ ■ 

6. CONCLUSIONS 

As we enter the Fermi era, it is worthwhile to make fi- 
nal predictions based on results from the EGRET CGRO 
era. Here we have adopted an approach making use of 
mean fluence ratios in the EGRET 100 MeV - 5 GeV and 
BATSE 20 kcV - 2 MeV bandpasses. Obviously this al- 
ready introduces significant K-corrections when calculat- 
ing the intrinsic properties of the GRB sources; a sample 
of GRBs with known redshifts is clearly preferred. 
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Fermi LAT (a and b) and AGILE (c and d) will detect for GRB 

Nevertheless, we fix the fluence ratios in these bands to 
obtain straightforward predictions for all-sky detection 
rates with Fermi, using on-axis effective areas and the 
BATSE 4B fluence distribution. Multiphcation by 1/5*'' 
gives the final de tection rate with the Fermi LAT tracker 
and calorimeter (jAtwood et al.l 12009*) 

For a 10% fluence ratio and a hard component with —2 
index, we estimate « 1 - 2 per month with more than 5 
photons with E > 100 MeV 7 rays in the LAT FoV, and 
1 or 2 per year with more than 100 7 rays with E > 100 
MeV. A bright GRB for spectral modeling, with more 
than « 100 photons with E > 100 MeV, is predicted 
once or twice per year. We note that with the use of —2 
fluence spectrum, the detection of a > 10 GeV photon 
from a GRB is considered very improbable, barring the 
existence of an additional component harder than —2. 

Furthermore, as the Fermi GBM and the Fermi LAT 
instruments begin to collect data from GRBs, we recom- 
mend to plot LAT fluence vs. GBM fluence and, more im- 
portantly, mean LAT GRB spectral index vs. LAT/GBM 
fluence ratios p to search for distinct classes of GRBs dis- 
tinguished by their high-energy properties. 

We thank David Band, Michael Briggs, Valerie Con- 
naughton, Brenda Dingus, Magda Gonzalez, Chryssa 
Kouveliotou, Julie McEnery and Nicola Omodei for 
many interesting discussions and useful suggestions, and 
especially to Nicola Omodei and Julie MeEnery for a 
really difficult and constructive report. We also thank 
the referee for comments, corrections, and useful sugges- 
tions. The work of CD is supported by the Office of 



GRB Predictions for the Fermi Telescope 



9 



Naval Research and NASA Fermi Science Investigation 
DPR-S-1563-Y, which also supported the research of TL 



at the Naval Research Laboratory. 



REFERENCES 



Atkins, W. B., et al. 2003, ApJ, 583, 824 
Atwood, R., et al. 2009, ApJ, to be submitted 
Band, D., et al. 1993, ApJ, 413, 281 
Band, D. L. 2002, ApJ, 578, 806 
Baring, M. G. 2006, ApJ, 650, 1004 

Belli, F. et al. 2007, Nuclear Instruments and Methods in Physics 

Research, A570, 276 
Briggs, M. S., et al. 1999, ApJ, 524, 82 

Casanova, S., Dingus, B. L., & Zhang, B. 2007, ApJ, 656, 306 

Catelli, J. R., Dingus, B. L., & Schneid, E. J. 1998, American 
Institute of Physics Conference Series, 428, 309 

Dermer, C. D., 2006, in The Multi-Messenger Approach to High- 
Energy Gamma-Ray Sources, eds. J. M. Paredes, O. Reimer, and 
D. F. Torres (New York: Springer), 127 

Dermer, C. D. 2007, ApJ, 659, 958 

Dermer, C. D., & Dingus, B. L. 2004, New Astr. Rev., 48, 537 
Dingus, B. L. 1995, Ap&SS, 231, 187 

Dingus, B. L., CateUi, J. R., & Schneid, E. J. 1998, American 

Institute of Physics Conference Series, 428, 349 
Gonzalez, M. M., Dingus, B. L., Kaneko, Y., Preece, R. D., Dermer, 

C. D., & Briggs, M. S. 2003, Nature, 424, 749 
Gonzalez, M. M., et al. 2009, ApJ, submitted 
Hurley, K., et al. 1994, Nature, 372, 652 

Kaneko, Y., Gonzalez, M. M., Preece, R. D., Dingus, B. L., & 

Briggs, M. S. 2008, ApJ, 677, 1168 
Kippen, R. M., et al. 2001, Gamma 2001: Gamma-Ray 

Astrophysics, 587, 801 
KouveUotou, C, Meegan, C. A., Fishman, G. J., Bhat, N. P., 

Briggs, M. S., Koshut, T. M., Paciesas, W. S., & Pendleton, 

G. N. 1993, ApJ, 413, LlOl 



Kouveliotou, C, et al. 1994, ApJ, 422, L59 
Kulkarni, S. R., et al. 1999, Nature, 398, 389 

Kwok, P. W., et al. 1993, American Institute of Physics Conference 

Series, 280, 855 
Le, T., & Dermer, C. D. 2007, ApJ, 661, 394 
Meegan, C. A., et al. 1996, ApJS, 106, 65 

Merck, M., et al. 1995, lAU CoUoq. 151: Flares and Flashes, 454, 
358 

Norris, J. P., & Gehrels, N. 2008, American Institute of Physics 

Conference Series, 1000, 280 
Omodei, N., & for the GLAST/LAT GRB Science Group 2006, 

astro-ph/0603762 
Paciesas, W. S., et al. 1999, ApJS, 122, 465 

Pittori, C, Tavani, M., & AGILE Team 2004, Bulletin of the 

American Astronomical Society, 36, 929 
Preece, R. D., Briggs, M. S., Mallozzi, R. S., Pendleton, G. N., 

Paciesas, W. S., & Band, D. L. 2000, ApJS, 126, 19 
Schneid, E. J., et al. 1992, A&A, 255, L13 
Schneid, E. J., et al. 1995, ApJ, 453, 95 
Sommer, M., et al. 1994, ApJ, 422, L63 
Sreekumar, P., et al. 1998, ApJ, 494, 523 

Stecker, F. W., Hunter, S. D., & Kniffen, D. A. 2008, Astroparticle 
Physics, 29, 25 

Strong, A. W., Moskalenko, I. V., & Reimer, O. 2004, ApJ, 613, 
956 

Thompson, D. J., et al. 1993, ApJS, 86, 629 

Thompson, D. J. 1986, NIM Physics Research A, 251, 390 

Wren, D. N., Bertsch, D. L., & Ritz, S. 2002, ApJ, 574, L47 



